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Abstract. We present results of Very Long Baseline Array (VLBA) observations of PKS 0528+134 at five frequencies (2.3, 
5.0, 8.4, 15.4, and 22.2 GHz). These quasi-simultaneous data enable us to study the spectral distribution of Very Long Baseline 
Interferometer (VLBI) components for the first time in this highly variable source, from which the central compact core is 
identified. Our observations indicate that there are two bendings for the jet motion at parsec scale. We provide an approximate 
spatial fit to the curved jet trajectory using the Steffen et al. 119951 helical model. We further investigate the proper motions 
of three jet components, which all show superluminal motion. At high frequencies (15.4 and 22.2 GHz) we detected a new 
component, which is estimated to be related to a radio burst peaking at about 2000. 
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1. Introduction 

As one of the brightest active galactic nuclei, PKS 0528+134 
has been intensively studied at the radio, X-ray, and gamma- 
ray bands. It is one of the furthest blazars detected by EGRET 
above 100 Mev (Mukherjee et al. 119971 with red-shift z=2.07 
(Hunter et al. 119931 . There are relatively less optical observa- 
tions compared with the observations at other bands because 
of its location close to the Galactic plane (galactic latitude 
b=-n°0'43") and thus high Galactic extinction. Whiting et 
al. d2"UU3l obtained B = 17.04 + 0.23, V = 16.74 ± 0.15 at 
2001.66 after the extinction corrections. But these are different 
from the previous values B=20.0 (Condon, Hicks & Jauncey 
IT9T71 and V=19.5 (Wall & Peacock fT9"83l . implying its optical 
variability. People have observed its variability on a timescale 
as short as a few days at gamma-ray (Mukherjee et al. I1996> . 
as short as two days in the optical band (Ghosh et al. 2000 1, 
and from several months to a few years in the radio band (e.g., 
Pohl et al. lT995l Peng et al. l200ll . There is a delay for the ra- 
dio burst from high frequencies to low frequencies. Peng et al. 
J200U showed in their light curves that the delay for the 1996 
burst is about 0.7 yr between 10.7 and 5.0 GHz. Mukherjee et 
al. 



served as a calibrator for VLBI observations (e.g. Shen et 
al.l2005l The first VLBI images of PKS 0528+134 (Chariot 
1990) show a one-sided extension at 2.3 GHz and a double 



19961 showed that the delay between the gamma-ray flare 
and the radio burst is a few months. 

Very Long Array (VLA) observations detected a diffuse 
emission 2" in extent and centered 1" north of the core with 
< 0.2% and 1.2% of the peak brightness at 4.9 and 1.5 GHz, 
respectively (Perlev [T9"g2l PKS 0528+134 is frequently ob- 



structure at 8.4 GHz, indicating a core-jet structure at PA. (po- 
sition angle) ~ -140°. Zhang et al. (1994) studied in detail 
the geometry, radiation mechanism, and physical environment 
of the emission region in the source through the X-ray and 
VLBI observations, and found a bent jet extending toward the 
northeast on a parsec scale. Krichbaum et al. ( 119951 hereafter 
K95) presented VLBI observations at 8, 22, 43 and 86 GHz, 
and also revealed a strongly bent one-sided core-jet structure 
with at least three superluminal components and two appar- 
ently stationary jet components. Their first 86 GHz VLBI ob- 
servations of this source revealed an unresolved core (<50/^as) 
and a new jet component (N2) whose ejection was believed to 
relate to an outburst in the mm-regime and a preceding gamma- 
ray flare observed in early 1993. Britzen et al. (1999, hereafter 
B99) studied in detail the morphology, kinetics, and the rela- 
tions between the morphological changes in the source struc- 
ture and the increased activity in the radio and gamma-ray 
bands through the eight-year 8 GHz geodetic VLBI observa- 
tions composed of 20 epochs. They found seven distinct jet 
components that all showed superluminal motions. The jet tra- 
jectory was bent twice, and different components had different 
paths. But these radio observations cannot give the spectral dis- 
tribution of VLBI components, because the source is strongly 
variable, while VLBI observations made at different frequen- 
cies were not simultaneous. 
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Quasi-simultaneous five-frequency VLB A observations of PKS 0528+134 



In Sect. 2, we present the first quasi-simultaneous VLBI 
mapping observations of PKS 0528+134 at five frequencies 
(2.3, 5.0, 8.4, 15.4, and 22.2 GHz). Then in Sect. 3 we discuss 
the structure of the jet, give the spectra of the jet components, 
identify the central core component based on its flat spectrum 
characteristic, study the proper motions of the jet components, 
estimate the ejection epochs of two components (a and n2), 
and describe the jet bending by means of the helical model, 
followed by the discussion and conclusion in Sects. 4 and 5, 
respectively. 

Throughout this paper, the radio spectral index a is defined 
as S v oc v~ a . By assuming Hq = 71 kms -1 Mpc~\ Qm = 0.27, 
and Q A = 0.73 (Spergel et al. l2003l . we have a scale for PKS 
0528+134 (z = 2.07) of 1 mas = 8.45 pc. 



2. Observations and data reduction 

The observations were performed with the NRAO 1 Very Long 
Baseline Array (VLBA) on August 20, 2001 (2001.64) at 
2.3/8.4 (dual frequency), 5.0, 15.4, and 22.2 GHz. Observations 
at different frequencies were interlaced, so data at different fre- 
quencies have similar (u,v) coverage. At 2.3/8.4 and 5.0 GHz, 
each scan lasted about 2.6 mins, with total observing time about 
36 mins each. At 15.4 and 22.2 GHz, scan lengths are about 25 
and 8 seconds with a total observing time of 60 and 38 mins, 
respectively. The data were recorded in 1-bit sampling VLBA 
format with a total bandwidth of 64 MHz (eight 8MHz IF chan- 
nels) at 5.0, 15.4, and 22.2 GHz, and the total bandwidth of 32 
MHz (four 8 MHz IF channels ) at the dual frequency 2.3/8.4 
GHz, per circular polarization at each station. 

The data correlation was done with the VLBA correla- 
tor in Socorro, New Mexico, USA. Post-correlation data was 
exported to the NRAO AIPS software (Schwab & Cotton 
119831 for fringe fitting, then to the Caltech DIFMAP package 
(Shepherd ll997> for hybrid mapping. A prior visibility ampli- 
tude calibration was done using the antenna gain and the sys- 
tem temperature measured at each station. During the hybrid 
mapping, the CLEAN and phase-only self-calibration were it- 
eratively used during the early processing, and the amplitude 
self-calibration was added for the later processing. The result- 
ing VLBA maps at five frequencies are shown in Fig. [0 The 
quantitative description of the source structure is determined by 
model fitting to the calibrated visibility data at each frequency, 
and the results listed in Tabled 

We used the Difwrap program (Lovell B000| to estimate 
the errors in Gaussian model fitting. The parameter of interest 
was stepped in small increments near the best-fit value. At each 
step this parameter was kept fixed, while some related parame- 
ters were allowed to relax during the model fitting. In addition, 
errors caused by self-calibration during the mapping processes 
had to be considered. These are typically 10% at 2.3, 5.0, 8.4, 
and 15.4 GHz, and about 15% at 22.2 GHz. The errors listed in 
Table[2were obtained by combining all these error budgets. 



1 The National Astronomy Observatory (NRAO) is operated by 
Associated Universities Inc., under cooperative agreement with the 
National Science Foundation. 



The single-dish flux density at the time close to our VLBA 
observations is 2.62±0.08, 2.40+0.08, and 2.52+0.04 Jy at 

5.0, 8.4, and 15.4 GHz, respectively, from the University of 
Michigan Radio Observatory, 2.20+0.03 Jy at 2.3 GHz from 
the NRAO Green Bank Interferometer, and 2.14 Jy at 22.2 
GHz from the Metsahovi Observatory. The total flux den- 
sity observed by these single dishes is 1.07+0.01, 1.06±0.03, 
1.06+0.04, 1.16+0.02, and 1.26±0.08 times higher than the in- 
tegrated flux density in the VLBA maps at 2.3, 5.0, 8.4, 15.4, 
and 22.2 GHz, respectively. We note that the relatively high 
ratio at 22.2 GHz was also seen at three epochs of VLBA ob- 
servations of PKS 0528+134 by Ojha et al. d2004t . These cor- 
rections to the absolute flux densities at different frequencies 
are applied in Sect. 3.2, where the spectral indexes of VLBI 
components are estimated. 

3. Results 

The main results are given in this section beginning with the jet 
components and an analysis of the jet morphology in Sect. 3.1. 
Then in Sect. 3.2, the core is identified through an analysis of 
the spectra of the jet components, and we also give the fitting 
to the spectrum of the core. In Sect. 3.3, we study the kinemat- 
ics of three jet components and discuss the ejections of the jet 
components. Last in Sect. 3.4, we try to use the helical model 
in S95 to describe the jet structure. 

3.1. Registrations of jet components and morphology 

After carefully analyzing the positions of each component at 
the five frequencies the registrations of the components are 
given in the first column of Tabled From Fig.^and Tabled 
we can see that component k in the south end is the most com- 
pact one and has the highest brightness temperature. B99 as- 
sumed it is the core, but we confirm it as the core from its spec- 
tral characteristics in the next subsection. The reason that nl 
component is not detected at frequencies higher than 2.3 GHz 
is probably due to its extended and weak structure. The sep- 
arations of component n2 from k at both 15.4 and 22.2 GHz 
are less than 0.3 mas. By considering that it is bigger than 0.3 
mas size at frequencies lower than 15.4 GHz, we think that k 
component at 2.3, 5.0, and 8.4 GHz in fact contains n2 com- 
ponent. Therefore the flux densities at 2.3, 5.0, and 8.4 GHz 
for k component actually include n2 emission. We treat n2 as 
a new component. Our VLBA maps at 15.4 and 22.2 GHz also 
indicate that there seems to be a new but quite weak compo- 
nent between components a and n2, about 0.65 mas from the 
core at PA. ~ 50°. However, the inclusion of this component 
affects the model fitting little, so we defer this to future obser- 
vations. We also note that the map at 5 GHz in Fig. 1 shows 
a very elongated component g. We tried to add an additional 
component in the model fitting. The fitted new component is 
very weak (~5 mJy) and has little effect on the fitting results. 
Thus, for simplicity we adopt component g only to represent 
the overall extended structure at 5 GHz (see Table 1). 

It can be seen from the maps that the overall morpholo- 
gies at each frequency are similar, they all extend in the north 
east direction with a moving tendency along the P.A.s;25 . It 
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Right Ascension (mas) Right Ascension (mas) 

Fig. 1. Naturally weighted maps of PKS 0528+134 at 2.3, 5.0, 8.4, 15.4, and 22.2 GHz. The fitting model components are shown 
with the corresponding names. The peak flux densities from low frequency to high frequency are 1.64, 1.97, 1.8, 1.65, and 1.26 
Jy/beam. The beam sizes and orientations from low frequency to high frequency are 6.77x4.04 mas 2 and -1.91°, 3.22x1.79 mas 2 
and -8.60°, 2.01x1.05 mas 2 and -7.78°, 1.1 1x0.54 mas 2 and -7.54°, 0.88x0.44 mas 2 and -6.1°. The off-source rms noise levels 
(cr) are 0.64, 0.51, 0.82, 1.0, and 2.1 mJy/beam. Contour levels are drawn at 3<x x(-l, 1, 2, 4, 8, 16, 32, 64, 128) at 22.2 GHz, and 
3cr x(-l, 1, 2, 4, 8, 16, 32, 64, 128, 256) at other frequencies. 



can also be seen from Table^that the jet's P.A. changes from 
about 50° at a separation from k component about 0.3 mas to 
40° at a separation about 1.4 mas and then to 25° at a separa- 
tion about 8 mas, clearly indicating the presence of jet bend- 
ing. K95 obtained maps at 8, 22, 43, and 86 GHz, which en- 
abled them to trace the jet down to sub-mas scale, and found a 
60-90° misalignment between the sub-mas and mas-structures. 
B99 found from their 8-year-long monitoring data of 8 GHz 
geodetic VLBI observations that the bending appears more pro- 
nounced near the core area (<1.5 mas) than farther away, and 
the jet appears to bend twice within the inner ~1.5 mas. These 
results are consistent with ours, except that our P.A. misalign- 
ment of 25° is less than that seen by K95. Combined with the 
data in Table [3] in B99, we can see that almost every compo- 
nent's P.A. decreases with the increasing separation from k un- 
til the P.A. reaches about 20° . We refer this to as the position an- 



gle regression. In Sect. 3.4 we will try to interpret the position 
angle regression using the helical model proposed by Steffen et 
al. (IT9951 hereafter S95). 

3.2. Identfication of the core and the spectra of other 
components 

The multi-frequency data can be used to study the spectra of 
the jet components. The spectra of the k component and other 
components are plotted in Fig.|3]and Fig. [3] respectively. Note 
that we have applied the absolute flux density correction when 
calculating the spectral indexes. There is a spectral turnover for 
the component k (Fig. Its spectral index between 15.4 and 
22.2 GHz is a — 0.30 + 0.08, much flatter than any other com- 
ponents (Fig. [3}- Component k is thus identified with the core. 
The identification of the core is important because it is thought 
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Table 1. The results of model-fitting. 



Component 


Flux(Jy) 


Separation(mas] 


i 0(deg) 


Major(mas) 


Axial ratio 


0(deg) 


T;,(k) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 








2.3 GHz, xi = 


1.02 








k 


0.532 ± 0.065 








0.52 + 0.10 


1.0 





1.39 x 10 12 


a 


1.114 ±0.11 


1.55 + 0.16 


32 + 3.3 


1.04 + 0.11 


1.0 





7.29 x 10" 


f 


0.318 + 0.034 


4.48 ± 0.56 


25 ± 2.6 


1.56 + 0.17 


1.0 





9.25 x 10 10 


g 


0.052 ± 0.006 


8.9 + 0.9 


27 ± 3.2 


3.40 ± 0.37 


1.0 





3.18 x 10 9 


nl 


0.022 ± 0.004 


19.7 + 2.0 


29 + 3.3 


8.80 + 0.91 


1.0 





2.01 x 10 s 








5.0 GHz, g = 


0.78 








k 


1.908 + 0.19 








0.54 ± 0.05 


0.50 + 0.01 


31 ± 3.1 


1.96 x 10 12 


a 


0.347 ± 0.036 


1.39 + 0.14 


41.2 + 4.1 


0.92 ± 0.09 


0.32 ± 0.02 


6 + 0.6 


1.92 x 10" 


f 


0.188 + 0.019 


3.65 ± 0.63 


25.0 ± 8.2 


2.7 ± 0.3 


0.57 + 0.01 


19 + 2 


6.78 x 10 9 


g 


0.027 ± 0.003 


9.18 + 0.93 


28 ± 2.8 


13.79+1.5 


0.16 + 0.02 


34 ± 3.4 


1.33 x 10 8 








8.4 GHz, xl = 


0.91 








k 


1.792 + 0.18 








0.33 ± 0.03 


0.43 ± 0.01 


47 ± 4.7 


2.03 x 10 12 


a 


0.390 ± 0.039 


1.22 + 0.12 


41.4 + 4.1 


1.70 + 0.17 


0.34 + 0.01 


28 ± 2.7 


2.11 x 10 10 


f 


0.069 ± 0.007 


4.2 ± 0.42 


24 ± 2.4 


1.78 + 0.18 


0.25 ± 0.07 


79+8.2 


4.62 x 10" 


g 


0.018 + 0.003 


8.4 ± 0.9 


22 ± 2.2 


4.21 ± 0.52 


1.0 





5.39 x 10 7 








15.4 GHz,x 2 v = 


: 1.57 








k 


1.405 + 0.14 








0.16 + 0.02 


1.0 





8.67 x 10" 


d2 


0.521 +0.052 


0.29 ± 0.03 


45 + 5 


0.26 ± 0.03 


1.0 





1.22 x 10" 


a 


0.179 + 0.018 


1.42 + 0.14 


42 + 4 


0.83 ± 0.08 


1.0 





4.10 x 10 9 


f 


0.064 ± 0.006 


3.59 + 0.36 


25 + 3 


1.75 + 0.18 


1.0 





3.30 x 10 8 








22.2 GHz, xl = 


: 1.83 








k 


1.156 + 0.17 








0.12 + 0.02 


1.0 





6.10x 10" 


n2 


0.347 ± 0.05 


0.28 ± 0.04 


50 + 8 


0.23 ± 0.03 


1.0 





4.98 x 10 10 


a 


0.104 + 0.016 


1.58 + 0.24 


39 + 6 


0.80 + 0.12 


1.0 





1.23 x 10 9 



Notes: Col. (1): component registration, col. (2): component's flux density in Jy, col. (3): component's separation from core 
component k in mas, col. (4): component's position angle relative to k from north to east in degree, col. (5): major axis 
of elliptical Gaussian component in mas, col. (6): ratio of the minor to major axis of elliptical Gaussian component, col. 
(7): position angle of the major axis of the elliptical Gaussian component from north to east in degree, col. (8): brightness 
temperature in K according to the formula given by Shen et al. (1997). 



to be related to the central engine in the AGN, and in general 
the jet is thought to be ejected from the core, powered by the 
central engine that is responsible for the activity of AGN. Other 
VLBI components (n2, a, f, and g) all have a typical, optically 
thin spectra (Fig.[3j, consistent with the core identification. The 
flux density of n2 at 2.3, 5.0, and 8.4 GHz, extrapolated from 
its spectral measurements at 15.4 and 22.2 GHz without con- 
sidering the spectral turnover, would be 4.76, 2.10, and 1.22 
Jy, respectively. At 2.3 GHz, this is much greater than the flux 
density of 0.53 Jy for component k, which is contradicted by 
the emission of n2 being included in the emission of k at 2.3 
GHz. Thus, the spectrum of n2 must be inverted at a frequency 
higher than 2.3 GHz. 

Furthermore, we simply try to fit the spectrum of the k 
component by using the synchrotron self-absorption (SSA) and 
free-free absorption (FFA) formula. According to Kameno et 
al. d2000T >. the SSA formula is 

S v = S v 25 [l - exp{-T s v- a - 2S \ (1) 

where So is the flux density in Jy at 1 GHz when the SSA opti- 
cal depth t s at 1 GHz » 1, and a is the optically thin spectral 
index. The FFA formula is 

S v = S()V~ a exp(-T f v~ 2A ), (2) 



where So is the flux density in Jy at 1 GHz extrapolated from 
the optically thin spectrum, and tj is FFA optical depth at 
1 GHz. Assuming a is known to be 0.30, we list the fit- 
ting results in Table [3] with the corresponding fitting curves 
also plotted in Fig. [2] The inverted spectrum can be fitted 
equally well by both absorption models. It can be seen in 
Fig. |2 that the spectral turnover happens at about 6.95 GHz 
with the corresponding flux density 1.93 Jy. Using the relation 
TTTT ~ 8( T 4,)5( 1 ^)§( 1 -^-)-5(l + z)s for the SSA model 

1 GHz M C v 1 Jy' y^mas' v ' 

(Kellermann & Pauliny 198 1 1, we can estimate a constant mag- 
netic field B = 0.35 mG in the core region. Here v m = 6.95 
GHz, S,„ = 1.93 Jy and 9 = 0.30 mas are used. Following 
Marscher (I1987> . we adopt 9 = y/9 a 9b as the angular diameter 
of the core, where 9 a and 9b are the major and minor axes of the 
Gaussian component. Since there is no VLBI size measurement 
at v„, = 6.95 GHz, the adopted 9 = 0.30 mas is an arithmetic 
average of the sizes at 5.0 and 8.4 GHz. The synchrotron cool- 
ing age is estimated to be at least 2.9 x 10 4 yrs using the expres- 
sion r 1/2 = 2.76 x 10 4 (yr) (^-^(^j-o s sin -i.s a (e . g ., 
Rybicki & Lightman ll979> . Here, B=0.35 mG, the turnover 
frequency in the source rest frame v^=21.33 GHz (after cor- 
recting the cosmological redshift effect), and the angle between 
the electron velocity and magnetic filed a- 1 . 
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Fig. 2. The fitted spectra of component k using the SSA (solid line) and FFA (dashed line) models. 
Table 2. The fitting parameters of Fig. 2. 



Model 


So(Jy) 


T 




1 

X 


SSA 


0.076 ± 0.006 


51 4 


5 


1.17 


FFA 


4.0 ± 0.2 


9.3 ± 


0.6 


1.78 



3.3. Proper motion and component ejection 

To find the relations between the seven components in B99 and 
the components in our paper, we extrapolate the proper mo- 
tions of the seven components in B99 to our observing time 
and find that the positions of a, f, and g components in B99 are 
near the positions of our a, f and g, so we think that our a, f, 
and g are the same as a, f, and g in B99. This is also the reason 
we used the same component registration. The counterparts to 
other components (b, c, d, and e) in B99 are not found in our 
VLBA observations. Figure @] shows the linear fit to the sepa- 
rations from the core of these VLBI jet components a, f, and g 
at 8.4 GHz. All the data before 2001.64 are from B99. Tabled 
lists the proper motion values of each component; due to the 
non-detections of b, c, d, and e in our observations, the proper 



motion values of components b, c, d, and e are from B99, but 
their apparent velocities are re-calculated using the new cos- 
mological constants (Ho - 71 kms _I Mpc~', Q.m = 0.27, and 
Qa = 0.73). It can be seen from Table that the proper mo- 
tion of the outmost component g (see Fig.Q becomes about 3 
times faster than those of other inner jet components a, b, c, d, 
e, and f, which are essentially the same about 10. 2-12. 7c within 
errors. 

Usually, the ejection of a new VLBI component can be re- 
lated to a flux outburst in the light curve (e.g. K95; B99), there- 
fore, we also relate the ejection of n2 to one burst. Assuming 
that there is no big difference in proper motion between n2 and 
other components, we adopt 0.13 mas/yr (see Table [5} as the 
proper motion for n2. Thus the ejection (zero-separation) epoch 
of n2 is roughly around 1999.5. This just falls into a burst du- 
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Fig. 3. The spectra of VLBI components in PKS 0528 + 1 34. Different lines represent the linear regressions performed on different 
VLBI components, i.e., n2, a, f, and g from top to bottom. 

Table 3. Combined with the data from B99, the proper motions are calculated assuming linear outward motion of the components. 



ID. 


a 




b 


c 




d 


e 


f 




T 


//(mas/yr) 

Papp 


0.13 ± 
11,0 1 


0.02 
1.7 


0.14 ±0.03 
11.8 ±2.5 


0.12 ± 
10.2 4 


0.01 
0.8 


0.13 ±0.01 
11.0 ± 0.8 


0.13 ±0.01 
11.0 ± 0.8 


0.15 ±0.02 
12.7 ± 1.7 


0.37: 

31.3 


I: 0.02 

± 1.7 



ration from 1999 to 2000.5 shown on the well-sampled light 
curves of PKS 0528+134 at 22 and 37 GHz (Terasranta et al. 
I2005> . suggesting that the ejection of n2 may be related to this 
radio burst, which peaks at about 2000. Likewise, the ejection 
epoch of component a is estimated to be about 1991.94+0.18 
from its proper motion 0. 13 mas/yr. This ejection epoch is con- 
sistent with the 1992.1+0.4 found by B99. Based on the ra- 
dio light curve at 8 GHz (Fig. in B99), B99 found that the 
ejection epoch of component a was near the beginning of the 
radio burst peaking around 1993.6, so they related the ejec- 
tion of component a to this radio burst, which is preceded by 
a gamma-ray burst peaking at about 1993.23 (Mukherjee et al. 
I1999i . But the light curves at 22 and 37 GHz (Terasranta et 



al. I2005i show another radio burst peaking at about 1992.32, 
and the flux-density variation at millimeter-bands shows a flux- 
density burst at about 1991.94 (K95), so we cannot exclude 
the possibility that the radio burst peaking at about 1992.32 or 
1991.94 may be related to the ejection of component a. The 
relations between the ejection of a new VLBI component and 
the flux-density bursts at radio and gamma-ray bands deserve 
more study to ascertain which flux-density burst is related to 
the ejection of the new VLBI component. 
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Fig. 4. The linear fitting to the separations of components a, f, and g from the core (component k). The data before 2001.64 are 
from B99, the ones at 2001.64 are ours at 8.4 GHz. 



3.4. Helical pattern of the jet structure 

VLBI observations have revealed a helical structures in some 
AGNs, e.g., 3C 345 (S95), 3C120 (Hardee et al. 120051 . and 
Mkn 501 (Villata & Raiteri [19991 . Two bendings of the jet in 
PKS 0528+134 on the mas scale indicate the possible presence 
of the helical trajectory in PKS 0528+134. The helical struc- 
ture may result from the orbital motion in a binary black-hole 
system, Newtonian precession, internal jet rotation, or global 
helical magnetic fields (e.g. Rieger l2~005> . S95 discussed the 
jet component's kinetics in detail through their helical models 
by using the conservation laws for kinetic energy, angular mo- 
mentum, momentum, and jet opening angle. 

Any three conservation quantities can give a jet trajectory, 
so there are four possible models. The first model assumes the 
conservations of momentum along the jet axis p z , the kinetic 
energy and the specific angular momentum L z , and gives a 
rotational hyperboloid with less than a quarter revolution. The 
second model assumes the conservations of p z , Ekin, and the 
opening angle if/, and gives a helical trajectory; but the helix is 



self-similar in space and time, and its helical amplitude is not 
dampened. The third model assumes the conservations of 
L z , and if/, and gives a helical trajectory that is not self-similar 
but has an asymptotic behavior; i.e., the pitch angle of the he- 
lix is the smallest at the origin, and it opens along the jet axis. 
The fourth model assumes the conservations of p z , L : , and iff 
and gives a similar trajectory to the one in the third model. The 
major difference between the third and fourth models lies in 
how the third model predicts a high level of component speeds 
as a component separates from the core. S95 applied the third 
model to fit the trajectories of C4 and C5 in 3C 345 simply be- 
cause this model reflects the basic character of the trajectories 
and motion observed in 3C 345. From Sect. 3.1 we know that 
the jet in PKS 0528+134 experiences two bends, and there is a 
position-angle regression phenomenon for the PA.s of the jet 
components, so obviously the first and second models are not 
applicable to the case of PKS 0528+134. From Tablef3]we can 
find that the apparent velocities of the jet components of PKS 
0528+134 remain on a high level as the components separate 
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from the core and appear to increase with distance from the 
core, so we reject the fourth model and adopt the third model 
to describe the bending trajectory of the jet in PKS 0528+134. 

The corresponding sets of equations for the third model in 
the cylindrical coordinates r,z, and <f>, are: 



r(t) = V(af + b) 2 + c 2 
r{t) - r 



z{t) 
<f>(t) 



tantft 



1 



sim// 



(arctan- 



at + b 



• arctan—) , 

c 



(3) 
(4) 

(5) 



where a = fisimf/, b 



^l-L 2 JP 2 



L z //3, p is the jet 



component's intrinsic velocity in units of the light speed, and 
ro,zo, and <po are the initial coordinates of the jet component 
(for simplicity, we set zo = 0). The best fit is not given here be- 
cause the currently available observational data points are not 
enough to constrain the evolution of the VLBI jet position with 
time, which requires 10 parameters to describe (S95). We show 
an approximate fit in Fig. [5] which is obtained through grad- 
ually adjusting the fitting parameters to minimize the differ- 
ence between the fitted trajectory and the observed data points. 
Some important fitting parameters are: the angle between the 
line of sight and the axis of cone is 3°, the intrinsic velocity 
P = 0.9995, the conical half opening angle = 0.25°, r = 0.1 1 
mas, and the initial angular velocity wo = -0.80°/yr. According 



to the relativistic beaming model, p app — 



PstnS 



, where p app is 



J "PP ~ l-/3cosa> "* l ^ L ^Happ 

the apparent velocity in units of the speed of light, and 6 is 
viewing angle to the moving direction of the jet. The allowable 
maximal value of 9 is 6 max = 2 arctan(;^-) (Pearson & Zensus 

Papp 

I1987> . so when p app -\3, 6> m£Ur =8.8°. Therefore, our choice of a 
viewing angle of 3° in Fig.[5]is reasonable. We can also esti- 



mate a Doppler factor from 5 



sinS 



When 6-3°, this gives 



<5 ~ 19, close to the value of 17 from 5 - 
is Lorentz factor) when using p 



T(\-pcos0) ^ -y 1-/3 2 

: 0.9995. We thus feel that the 
fit shown in Fig.|5]is physically reasonable. 



Furthermore, this model can be used to explain the 
position-angle regression discussed in Sect. 3.1. Since L z = rfi^ 
is conservative (/?„ is the velocity of the jet component rotating 
around the conical axis), p ia is decreasing with the increasing 
r when the jet component is moving out along the conical sur- 
face. By differentiating Eq. (3), one can obtain the radial ve- 



locity component p r - 



a(at+b) 



:, which is perpendicular to p. 



-\J(at+b) 2 

and p^ and which increases with the time. Similarly, from Eq. 
(4), one can get p : = jjj^, which increases much faster than p r 
as tant/f «: 1; i.e., p u is mainly converted to p z . There is thus 
a reacceleration process of the velocity component p : moving 
along the jet axis, and finally the motion of the jet component 
is dominated by p z . Thus, the jet component moves just about 
along the conical axis, and the PA. barely varies when the jet 
component moves into the outer region of the jet; i.e., the mov- 
ing direction on the map of the jet component regresses to the 
projected direction of the cone axis on the sky. 



4. Discussion 

Our quasi-simultaneous VLBI observations at five frequencies 
allow us to study the spectral distribution of PKS 0528+134 
for the first time. From the light curves at 22 and 37 GHz 
(Terasranta et al. |2005 1, we can see that our observing time 
2001.64 PKS 0528+134 is in the quiescent state, so the spec- 
tral distribution represents the spectral distribution of PKS 
0528+134 in its quiescent state. The identification of the core 
in VLBI analysis is important because the registration, tracing, 
and kinematic analysis of the jet components are all based on 
such a core identification. B99 assumed that the brightest and 
most compact component at the southern end of the jet was 
the core, based on how all other VLBI components seemed to 
move away from it, and how it showed a pronounced flux den- 
sity variation and an inverted spectrum (Pohl et al. 1995). 

Our spectral analysis indicates that component k has a flat 
spectral index, a = 0.3 + 0.08, with the spectrum inverted at 
6.95 GHz, consistent with the k component being the com- 
pact core. According to the SSA theory, the separation of an 
optically-thin jet component to the core increases with the in- 
creasing frequency (Lobanov; 1998 1. However, due to large un- 
certainties in the position determination in our observations, we 
cannot tell whether there is such a core-position offset in PKS 
0528+134. Phase-referenced VLBI observations can achieve a 
high positional precision of tens of microarcsecond (Ros 2003 1 
and can thus be used to study the absorption mechanism in 
the core region of PKS 0528+134 through investigating the 
frequency-dependent position offset. 

By tracing the components in B99, we find our a, f, and 
g components are identical to the a, f, and g components in 
B99. All three components show high apparent velocities, es- 
pecially g up to 31.3c. The ejection times of component a and 
the newly detected n2 component can be estimated from the ex- 
trapolation of the proper motions to their zero separations from 
the core. The ejection of n2 at about 1999.5 could be related 
to the radio burst peaking at about 2000, but we cannot be sure 
which flux-density burst is related to the ejection of component 
a. We can see from the light curves at 22 and 37 GHz given 



by Terasranta et al. |20O5 1 that there are six significant flux- 
density outbursts before our observation at epoch 2001.64, but 
only two outbursts could be related to the ejections of VLBI 
components. The reason that not every flux density outburst 
corresponds to an ejection of a VLBI component could be that 
the related VLBI components are too weak or have reached its 
radiation loss time before our observations. Also it is not so 
clear whether or not every outburst can be related to the ejec- 
tion of a VLBI component. 

According to the morphology of PKS 0528+134, we as- 
sume for the first time that it has a helical structure. New 
VLBA observations at 15 GHz (Lister & Homan I20051 have 
also confirmed the presence of a bending jet structure in PKS 
0528+134. We used the third helical model in S95 to describe 
the jet trajectory in PKS 0528+134. S95 fitted this helical 
model to the trajectories of C4 and C5 in 3C 345. In 3C 345 
there is a component, 18 mas from the core, whose PA. is dif- 
ferent from the PA.s of its inner components by 20° (Unwin 
& Wehrle [T992> . which is contrary to the general prediction of 
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Fig. 5. We use the third model of Steffen et al.'s (1995) helical models to approximately fit the trajectory of the jet components 
obtained from our quasi-simultaneous five-frequency VLBA maps. The locations of the components relative to the core (compo- 
nent k in Fig. 1) are the average values of the locations measured at different frequencies. The standard deviation of the locations 
of each component relative to the average is adopted as the error bar. 



this helical model that the jet component will move close to 
the jet axis little by little; i.e., there is a position-angle regres- 
sion. Because of this, S95 only restricted their analysis to the 
inner 3 mas from the core. Our VLBA observations at 2.3 GHz 
(Fig-D show that there is no obvious P.A. change for the exter- 
nal jet components of PKS 0528 + 134. We then can apply this 
helical model to the jet of PKS 0528+134 until 20 mas from 
the core since there is no abrupt P.A. change for the external 
jet components as seen in 3C 345. Some core-dominated radio 
sources even show 90° misalignment between the PA.s of the 
VLBI and VLA scale jets (Conway & Murphy [T9"9"3T l. which is 
difficult to explain by the helical model we adopted. Conway 
& Murphy d!993> used their low-pitch helical model to ex- 
plain this phenomenon, but this model is difficult to explain 



the several bends at the VLBI scale in radio sources like PKS 
0528+134 and 3C 345. We think that the tremendous change 
of the P.A. on a large scale may be due to the interactions be- 
tween the jet and ambient medium of the host galaxy or another 
external force performed on the jet, so the angular momentum 
conservation required by the helical model we used in Sect. 3.4 
is violated. Additional dynamical analysis would be required 
in order to study the large-scale jet structure. Except that many 
extragalactic radio sources have the helical jets, recent VLBI 
observations show that some microquasars in our Galaxy also 
have these helical patterns. For example, both GRO J1655-40 
(Hjell ming & Rupen fT995l and SS433 (Blundell & Bowler 
I2004t exhibit the helical jets. 
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We know from the analysis of the helical model that the 
velocity of jet component is dominated by the velocity com- 
ponent along the cone axis, i.e., fi z . According to the helical 
jet model, {S~ should first increase rapidly and then approach its 
theoretical value gradually if the angle between the line of sight 
and the z-axis is constant. From Fig.|4]we find it impossible 
that the apparent velocities of components a and f are reaccel- 
erated to that of component g by this helical model. A more 
realistic situation could be that the velocity component {S z s of 
these jet components have been accelerated shortly after their 
emergence, and that the initial ejected velocity of g component 
is bigger than others'. Obviously, this scenario requires a very 
good time coverage of high-resolution VLBI observations to 
test it, especially during the period when the jet component is 
initially ejected. 

5. Conclusions 

The five high-resolution VLBI maps made at frequencies 2.3, 
5.0, 8.4, 15.4, and 22.2 GHz show that PKS 0528+134 has a 
jet structure extending in the north east direction, and the jet 
bending is clearly seen at about 0.3 and 1.4 mas from the core. 
Our registration of VLBI components a, f, and g is consistent 
with the one in B99. We also detect a new VLBI jet component 
n2, which is likely to be ejected during the radio burst at about 
1999.5. We estimate the ejection epoch of component a to be 
1991.94 by the linear back-extrapolation of its proper motion, 
but cannot determine which flux-density burst is truly related to 
the ejection of component a. No counter-jet is seen in our maps, 
but its existence cannot be ruled out. We also find that the po- 
sition angles of the jet components in PKS 0528+134 are large 
when they are initially ejected from the core, but gradually de- 
crease with the increasing separation from the core until about 
25°, which is referred to as the position angle regression. 

This first quasi-simultaneous VLBA observations at five 
frequencies provide a good condition to calculate the spectra 
of the jet components without taking the variations of the flux 
densities with the time into account. We identify k component 
as the core mainly based on its flat spectral index of 0.30 + 0.08 
and an inverted spectrum peaking at about 6.95 GHz. But we 
cannot determine which absorption mechanism (SSA, FFA, or 
others) is mainly responsible for the observed inverted spec- 
trum; future low-frequency observation with sufficient spatial 
resolution may be useful for distinguishing these possible mod- 
els. The brightness temperature Tj, of each component is given 
in Table [fl the brightness temperature of core component k 
does not exceed the equipartition brightness temperature cutoff 
of ~ 10 11 K (Readhead ll994> if we take the Doppler boosting 
effect into account. The presence of fast superluminal motions 
indicates that the relativistic effects are important in this source. 

We attempt to describe the trajectory of the jet in PKS 
0528+134 using the helical jet model proposed by S95. The 
parameters of kinetics obtain from this helical model are con- 
sistent with ones from observations. We find from this helical 
model that there is a reacceleration process of the velocity com- 
ponent/?; along the jet axis and that the jet finally moves nearly 
along the jet axis, which can be used to explain the position- 
angle regression. Since the reacceleration process in this helical 



model cannot explain this big difference in the jet's velocities, 
we suggest further that a much higher apparent velocity of g 
component than its inner components' may be due to the rela- 
tively higher initial velocity of g component than its inner com- 
ponents. To understand the reacceleraton process well in this 
helical model needs a good time coverage of high-resolution 
VLBI observation just as the jet component ejects. 
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